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Abstract. Perturbative QCD in the high-energy hmit describes the evolution of scattering ampH- 
tudes with increasing energy towards and into the so-called saturation regime. Comparisons of the 
predictions with experimental data for a number of observables led to significant progress and un- 
derstanding. We discuss the case of particle-production cross-sections measured at HERA and argue 
that these measurements have the potential to provide evidence for the saturation regime of QCD. 

INTRODUCTION 

In the Regge limit of perturbative QCD, i.e. when the centre-of-mass energy in a col- 
lision is much bigger than the fixed hard scale of the problem, parton densities inside 
the colliding particules grow with increasing energy, leading to the growth of scattering 
amplitudes. When the parton density becomes too large and the scattering amplitudes 
approach the unitarity limit, one enters in a regime called saturation ||l|]. 

The transition to the saturation regime is characterized by the so-called saturation 
scale which is an intrinsic hard scale of the problem. Contributions to the scattering 
amplitudes which are neglected as higher twist in the Bjorken limit of perturbative QCD 
become important in the saturation regime: leading-twist gluon distributions are no more 
sufficient to describe scattering at high energies. 

A consistent approach is to express physical observables in terms of the leading 
terms in an expansion with respect to the inverse of the center-of-mass energy: this is 
called the eikonal approximation. This formalism is well-suited because, as the energy 
increases, density effects and non-linearities that lead to saturation and unitarization of 
the scattering amplitudes can be taken into account. 

In the following, after introducing the eikonal approximation, we discuss the case 
of particle-production cross-sections measured at HERA. We concentrate on the phe- 
nomenology for diffractive observables: we review their model descriptions and investi- 
gate their potential to provide evidence for parton saturation. 

HIGH-ENERGY SCATTERING AND SATURATION 

Let us start with the eikonal approximation for quarks and gluons scattering at high 
energies. When a system of partons propagating at nearly the speed of light passes 
through a target and interacts with its gluon fields, the dominant couplings are eikonal: 
the partons have frozen transverse coordinates and the gluon fields of the target do not 
vary during the interaction. This is justified since the time of propagation through the 




FIGURE 1. Eikonal propagation of a quark with transverse position x through a target evolved at 
rapidity Y. The eikonal phase Wf{x), see formula resums all number of gluon exchanges. 



target is much shorter than the natural time scale on which the target fields vary. The 
effect of the interaction with the target is that the partonic components of the incident 
wavefunction pick up eikonal phases: if \{a,x)) (resp. \{a,x))) is the wavefunction of 
an incoming quark of color index a E [l,A/c] (resp. gluon of color index aE[l,N^ — l]) 
and transverse position x (the irrelevant degrees of freedom like spins or polarizations 
are not explicitly mentioned), then the action of the 5— matrix is (see for instance [2]): 

5|(a,x))®|r) = £[W'^(x)]„„,|(a',x))®|0, S\{a,x))(»\t) = Y,Wf{x)\{b,x))®\t), 

a' b 

(1) 

where \t) denotes the initial state of the target. The phase shifts due to the interaction are 
the color matrices Wf and Wa, the eikonal Wilson lines in the fundamental and adjoint 
representations respectively, corresponding to propagating quarks and gluons. They are 
given by 

WfA^) =Pexp{ig,J dz+T^^AA"_{x,z+)} (2) 

with A_ the gauge field of the target and Tp^ the generators of SU (Nc) in the fundamen- 
tal (F) or adjoint (A) representations. We use the light-cone gauge = and P denotes 
an ordering in the light-cone variable z+ along which the incoming partons are propa- 
gating. As displayed in Fig. 1, all number of gluons exchanges are included in ([T]), which 
shows the leading term of an expansion with respect to the inverse of the center-of-mass 
energy squared 5 ~ where Y is the total rapidity. 

For an incoming state |^,>,), the outgoing state \^out)=S\^in)®\t) emerging from the 
eikonal interaction is obtained by the action of the 5— matrix on the partonic components 
of \^in) as indicated by formula ([T]). The outgoing wavefunction \^out) is therefore a 
function of the Wilson lines Q. When calculating physical observables from \^out), 
one obtains objects that are target averages of traces of Wilson lines (the traces come 
from the color summation). For instance, the simplest of these objects is 

Tgg{x,x';Y) = 1 - ^ (tv(w^{x')Wf{x)') )^ , (3) 

namely the dipole scattering amplitude (jc, x' : positions of the quark and antiquark) 
off a target evolved at rapidity Y. The target average has been denoted { .)y and contains 
the Y dependence. The amplitude Q enters for instance in the DIS total cross-section 



(see next section). More generally, observables are functions of © or more complicated 
amplitudes. Let us introduce another one of them, which we shall need later: 

T^f{x,x';y,y'-Y) = 1 - ^ {Tr{w^{x')WF{x))Tr{w^iy')Wfiy)) )^ . (4) 

This is the scattering amplitude for a set of two dipoles (x,x') and Cv,/). The amplitudes 
© and @ take values between (transparency) and the black-disk (saturation) limit 1. 

To actually compute these amplitudes, one has to evaluate the averages { ■ )y vvhich 
amounts to calculating averages of Wilson lines: = / ^A/[A]i7y[^] where the 

target wavefunction Uy [A] represents the probability to find a given field configuration 
inside the target evolved at rapidity Y. The information contained in the target averages 
is mainly non-perturbative but the evolution towards higher rapidities dUY[A]/dY can 
be computed perturbatively, at least in the leading-logarithmic approximation. Several 
equations have been established with different degrees of approximations, we shall not 
discuss them here and the reader can refer to |3] for more details. Let us only mention 
the Balitsky-Kovchegov saturation equation (BK) ^ which is a closed equation for T^q 
obtained in a mean-field approximation. We shall refer to the BK equation later on when 
we link observables to the dipole amplitudes © and @ and discuss phenomenology. 

SATURATION PHENOMENOLOGY AT HERA 

In deep inelastic scattering (DIS), a photon of virtuality collides with a proton. In an 
appropriate "dipole" frame, the virtual photon undergoes the hadronic interaction via a 
fluctuation into a dipole (see Fig. 2, left); the dipole then interacts with the target proton 
and one has the following factorization 

(yDis{Q\Y) = jd^rxt/{\r\,Q^)2jd^bT,q{b-^,b+'-;Y^ (5) 

which relates the DIS total cross-section 00/5 to the qq— dipole amplitude Tqq. The 
function V^(r, 2^) = J dz |(^^(r,z; 2^) |^ is obtained from the well-known wavefunction 
^^{f,z',Q^) which describes the splitting of the photon onto a dipole of transverse size 
r and with the antiquark carrying a fraction of photon longitudinal momentum z- Note 
that in this case, not all the information on T^q is relevant as the impact parameter b is 
integrated out: only the cross-section Oqq{r,Y) =2 J d^b Tqq{b—^,b+^\Y^ is needed. 

Measurements of Odis at HERA have had a great impact on saturation phenomenol- 
ogy, especially the discovery of geometric scaling [5]: the fact that OoisiQ^i Y) is a func- 
tion of the single variable / QI{Y) with the saturation scale 2^(7) ~exp {XY) and A ~ 
0.28. Indeed this has a natural explanation in terms of traveling-wave solutions [6] of the 
BK equation Tqq{r,b = 0,Y) = Tqq{rQs{Y)). However, this result is obtained neglecting 
the impact parameter dependence of Tqq and considering Oqq {r,Y)=Spx Tqq (r, Z? = 0, 7 ) 
where Sp is the transverse area of the proton fitted to the data. 

In order to understand better and study more consistently the impact parameter de- 
pendence of Tqq, the authors of [7] have looked at diffractive vector-meson production 
(see Fig. 2, center). In diffractive deep inelastic scattering, the proton gets out of the col- 
lision intact and there is a rapidity gap between that proton and the final state X. When 



FIGURE 2. Three processes measured in virtual photon-proton colHsions at HERA: DIS total cross- 
section (left), diffractive vector-meson production (center), and diffractive photon dissociation (right). 



the final state is a vector meson, the momentum transfer A dependence of the cross- 
section is related to the impact parameter dependence of the dipole amplitude. Indeed 
the cross-section reads = — A^) 



do 1 

dt ~ An 



I d\ ^>{\rlQ\Ml) I d^b e*-^ T,^(^b-^,b + ^;Y 



(6) 



where the function ^(r,^^, My) = J dz ^'^{r,z;Q^)^^ {r,z;Mv) is obtained from both the 
photon wavefunction 0^ and the final-state meson (whose mass has been denoted My) 
wavefunction 0^. By analysing data on p— meson production at fixed and 7 ~7, the 
authors extracted the dipole 5— matrix Sqq{r,b;Y) = l — Tqq{b—^,b-\-^;Y^ as a function 
of b for a fixed size tq with ''g ~ 4/ (Q^+My). Three different sets of data at different 

have been used. Their results are shown on the left plot of Fig. 3; the shaded area on 
the left is an uncontrolled region due to the lack of large—? data. The plot shows that the 
b dependence cannot be neglected and that Tqqijr^ 1 GeV~\Z?~0;F~7) ~0.4. This 
significant value of Tqq indicates that HERA could be entering the saturation regime. 

As the importance of the impact parameter had been pointed out, a phenomenological 
model for the dipole amplitude Tqq with an impact parameter profile was proposed in 
fS]. With that saturation parametrization, the authors could well reproduce the data 
for diffractive J /m production at HERA: the t spectrum as well the the and Y 
dependences of the total cross-section o>/ip. Their results are displayed on the center 
and right plots of Fig. 3 where one can see the good agreement of the model with the 
data. A successful description of the same data using numerical simulations of the BK 
equation was also given in confirming the compatibility of saturation predictions. 

However in all the model descriptions of t spectra, the impact parameter depen- 
dence was introduced by hand as one had not extracted any information on the b de- 
pencence of Tqq from saturation equations. That moderated the impact of the results 
mentioned above. Interestingly, it was recently O] pointed out that important in- 
formation can be obtained from the BK equation when looking at the A dependence of 
7'(r,A;y) = / d^b e'*'^ Tqq(^b—^,b+^;Yy For instance, the geometric scaling property 
was extended: f{r,A;Y) = f{rQs{A,Y)) with 2^(A, 7) ^A^exp (AY). Parametrizing the 
dipole amplitude with the momentum transfer instead of the impact parameter opens a 
new approach to analyse the data. An experimental confirmation of geometric scaling at 
non-zero momentum transfer would represent a significant success for saturation. 



y'p -> J/* p 




FIGURE 3. Left plot: the q'^— dipole 5— matrix extracted from the diffractive p —meson production data 
in 0]; this shows the impact parameter dependence for three different dipole sizes rg and a rapidity 7 ~7. 
Center and right plots: diffractive 7/^ production at HERA; da/dt as a function of f (center), CT/yij/ as a 
function of W ^ exp(F /2) (top right) and (bottom right); comparison with the model of ||8|] is shown. 



Let us finally consider high-mass diffraction. If the final state diffractive mass Mx is 
much bigger than Q, the dominant configurations to the final state come from the qqg 
component of the photon wavefunction (see Fig.2, right) or from higher Fock states, 
i.e. from the photon dissociation. By contrast, if Mx ^ 2, the dominant configurations 
come from the qq component as it was the case for vector-meson production. Let us then 
consider the kinematical regime where P = Q^/ (Q^+M^) <^ 1 and investigate the qqg 
component. The right plot of Fig.2 represents the diffractive production of a gluon with 
transverse momentum k and rapidity log(l//3) in the collision of the photon with the 
target proton. Provided ^ is a hard scale, the gluon momentum spectrum is given by [|12I1 



Mxdo ttsN- 



(7) 



(fkdMx Itz'^Cf 

where At] = 7 — log( 1 //3 ) is the rapidity gap. The two-dimensional vector A is given by 



A{k,x,x';An)=J^e-''-' (r^f {x,z;z,x';An) -T,^{x,x'-Ari)) ■ 

(8) 

(2) 

Interestingly enough, independently of the form of the dipole amplitudes Tqq and T^y , 

the behavior of the observable k^ do / d^kdMx as a function of the gluon transverse 
momentum k is the following lUsIl : it rises as k^ for small values of k and falls as l/k^ 
for large values of k. A maximum occurs for a value which is related to the inverse of 
the typical size for which the T— matrices approach one; in other words, the maximum 
ko reflects the scale at which unitarity sets in. If the energy is large enough so that the 
saturation scale is hard, unitarity will come as a consequence of parton saturation and 
ko ~ Qs- If not the case, unitarity will be rather due to non-perturbative physics. 




FIGURE 4. The diffractive gluon production cross-section k^Mxda /dk^dMx- Left plot: as a function 
of the rescaled gluon transverse momentum k/Qs for two extreme values of equal to 0. 1 and 100 GeV^ 
and four values of the saturation scale = 0.5, 1,2,3 GeV. Right plot: as a function of the jet transverse 
momentum k and in the HERA energy range for 2^ = 0.1 GeV^ and different values of diffractive mass 
Mx and energy W; full lines: only the light quarks are included in i//, dashed lines: charm is also included. 



In the saturation case, the model of for the dipole amplitudes allows to plot the 
whole k spectrum O. This is shown on the left plot of Fig.4 and one can indeed see 
that the spectrum features a maximum peaked around — 1 -4 Qs independently of 2^ 
and Qs. Measuring this cross-section at HERA would offer a unique opportunity to test if 
saturation plays a role in diffraction at the present energies. On the right plot of Fig.4, the 
cross-section is plotted in the HERA energy range for different values of Mx and total 
energy corresponding to different values of Qs- The saturation scale is the one 

extracted |[l4ll from the F2 data. As expected for realitic jet transverse momenta, k>Qs 
and the data would lie on the perturbative side of the bump. There is a big difference in 
the rise towards the bump between the lowest (top right) and highest (bottom left) 
bins. A confirmation of this behavior would certainly favor the saturation scenario. 



REFERENCES 

1. L. V. Gribov, E. M. Levin and M. G. Ryskin, Phys. Rep. 100 (1983) 1. 

2. A. Hebecker, Phys. Kept. 331 (2000) 1; A. Kovner and U. Wiedemann, Phys. Rev. D64 (2001) 1 14002. 

3. L. McLerran, these proceedings; R. Peschanski, these proceedings. 

4. I. BaHtsky, Nucl. Phys. B463 (1996) 99; Y. V. Kovchegov, Phys. Rev. D60 (1999) 034008. 

5. A. M. Stasto, K. Golec-Biernat and J. Kwiecinski, Phys. Rev Lett. 86 (2001) 596. 

6. S. Munier and R. Peschanski, Phys. Rev Lett. 91 (2003) 232001; P/zy,?. Rev D69 (2004) 034008. 

7. S. Munier, A. M. Stasto and A. H. Mueller, Nucl. Phys. B603 (2001) 427. 

8. H. Kowalski and D. Teaney, Phys. Rev D68 (2003) 1 14005. 

9. E. Gotsman, E. Levin, M. LubHnsky, U. Maor and E. Naftali, Acta Phys. Polon. B34 (2003) 3255. 

10. C. Marquet, R. Peschanski and G. Soyez, Nucl. Phys. A756 (2005) 399. 

11. C. Marquet and G. Soyez, Nucl. Phys. A760 (2005) 208. 

12. C. Marquet, Nucl. Phys. B705 (2005) 319; Nucl. Phys. A755 (2005) 603c. 

13. K. Golec-Biernat and C. Marquet, Phys. Rev D71 (2005) 1 14005; hep-ph/0509034 

14. K. Golec-Biernat and M. Wiisthoff, Phys. Rev D59 (1999) 014017; Phys. Rev D60 (1999) 1 14023. 



